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Abstract. Waves play a crucial role in diagnosing the plasma 
properties of various structures in the solar corona and coronal 
heating. Slow magneto-acoustic (MA) waves are one of the im- 
portant magnetohydrodynamic waves. In past decades, numer- 
ous slow MA waves were detected above the active regions and 
coronal holes, but rarely found elsewhere. Here, we investigate a 
'tornado'-likc structure consisting of quasi-periodic streaks within 
a dark cavity at about 40-110 Mm above the quiet-Sun region on 
2011 September 25. Our analysis reveals that these streaks are ac- 
tually slow MA wave trains. The properties of these wave trains, 
including the phase speed, compression ratio, kinetic energy den- 
sity, etc., are similar to those of the reported slow MA waves, 
except that the period of these waves is about 50 s, much shorter 
than the typical reported values (3-5 minutes). 

1 Introduction 

Magnetohydrodynamic (MHD) waves are important phe- 
nomena in solar corona and catch broad interests. They are 
believed to be the energy carriers for coronal heating [e.g. , 
Erdel vi and Fedun\ , [200 7j; \Jess et oil 120091 : \McIntosh et all 

diagnose plasma 



to 



201l||, and also be a powerful tool 
prope rties and physi c al processes in t he sol a r atmosphere 
[e.g., \ Robert! l200Cl : I Baneriee et all l2007t I Marsh et all 
[2003|. Slow magneto-acoustic (MA) waves are one of the 
eigenmodes of MHD waves, and were frequently observed 



above the active regions, e.g., in coro nal loops, and above 



coronal holes, e.g., in plum es [e. g.. \Ofrnan et all Il997 



1999 



199E ; \DeForest and Gurman L 19981 : B^i^fmiajis^idClette 



Berahrnans et ai , 200ll: \Scliriive^^t~cd[ , 



Ni^ktxn^^^ al^ 19991: I De Moortel et 



Baneriee et all 120071 : \Mariska and Mualachl l20l"d] 



The 

propagation speed of these waves is usually around 100 
km s~ x with a period of about 3 min above active regions 
or of about 10 — 15 min above coronal holes. Due to 
the different magnetic field configurations, they could be 
found over a wide range of altitudes from 7 Mm to about 
several hundred Meg a meters above coronal h oles [e.g., 
\DeForest and Gurman] , Il998l : \Ofman et all . Il997t [, but usu- 
ally below 20 Mm in coronal l oops above active regions [e.g., 
\Berahmans and ClettA , 1 19991 : \Pe MooHel et all l2002bl la| . 



In past decades, slow MA waves above quiet Sun re- 
gions were rarely reported. It was probably due to the 
limitation of instruments. Before t he era of the Sola r 
Dynamic Observatory (SDO) mission \Pesnell efaH l2012f | . 
imaging instruments have provided either full-disk observa- 
tions at relatively low spatial resolution and low cadence 
(e.g., SOHO/EIT referring to \Pelaboudiniere et q/.l Il995l . 
STEREO/EUVI, I Wuelser et all 12004 ). or high resolution 
and high cadence but limite d field of vie w (FOV) (e.g., 
TRA CE, \Handv et al\ Il999l , Hinode/EIS , \Culhane et al\ 
Il997f ). The instrument AIA \hemen et all |2012| | on board 
SDO provides so far the best combination of resolution, ca- 
dence, FOV and EUV wavelength coverage. In particular, its 
full-disk and round-the-clock coverage allows one to launch 
a comprehensive search for a phenomenon of interest. The 
only shortcoming of AIA data is the lack of Doppler in- 
formation. Although this shortcoming brings difficulty in 
the identification of waves, it could be compensated by the 
analysis of multi-wavelength observations. In this paper, we 
report an observation of multiple slow MA wave trains in 
a dark cavity above a quiet-Sun region, which may advance 
our understanding on the MA waves in the solar atmosphere. 
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2 Overview 

The event of interest took place on 2011 September 25 
above the southwest limb within a dark cavity above the 
quiet Sun region (as shown in Fig. [TJa)). The overall cavity 
structure stayed stable all the time, but the prominence-like 
material under the cavity was perturbed at the beginning 
of September 25 and finally settled down on the next day. 
The process from 08:00 to 13:00 UT on September 25 can be 
seen in the o nline movie M l. Among various phenomena in 
this process, \Li et al\ [20121 ] reported a so-called 'huge solar 
tornado'. 

In this study, we focus on the nature of the repeated quasi- 
periodic (QP) streaks starting from about 9:00 UT. The 
most clear group of QP streaks (labeled as SI in Fig[lld)) 
appeared around 10:10 UT in the 171 A channel, whose 
peak formation temperature is about 0.63 MK. Several 171 
A images taken before, during and after this interval are dis- 
played in Figure [T] Before the appearance of the QP streaks, 
there were apparently two plasma tubes (labeled as 1 and 2) 
standing right above the limb in the dark cavity (Fig-D^b)). 
Initially, the two tubes were slightly tangled as if they were 
crossed each other in projection. Thereafter, those two tubes 
were disturbed by an Ml. 5 flare in the nearby active region 
NOAA 11303 (Fig.[TJa)) and ended up 'apparently' merging 
together after 9:20 UT. 

Meanwhile, continuous upward flows along Tube 1 can 
be noticed to originate from the tube bottom starting from 
around 08:20 UT. The upward flows seemingly transformed 
into helical-like flows with increasing height, and near the 
top of Tube 1 the helical-like flows quickly expanded and 
developed into QP streaks, which intermittently appeared 
with various intensity and broadness in the next 4 hours. 
Figure [He) snows another group of streaks (labeled as S2). 
These streaks propagated along the tube once they ap- 
peared. 

3 Helical magnetic structure, intermittent 
mass flows or wave trains? 

The nature of the QP streaks is of particular interest. In- 
tuitively, one may think that these QP streaks reflected the 
helical motion of plasmas in a magnetic flux rope. If this 
conjecture is true, the degree of the twist of the magnetic 
field in the flux rope could be inferred from the number of 
streaks. In the 171 A image at 10:10 UT, for example, at 
least 10 streaks were displayed at the upper portion of the 
tube, and it implies that the magnetic field lines twisted 
more than 10 rounds or 207T in radian. Previous works have 
suggested that a flux rope with force-free field tends to be 
kink unstable when the twist exceeds ab out 3.57T, and the 
grow t h rate is on the order of minutes \Hood and PriesA 
I1979I : \Tordk and Klierrl 120031 : \T6rdk et all 120041 ] . Some re- 
cent studies further suggested that, for a flux rope with 
a larger aspect ratio or axial mass flow, the threshold 
may increase up to 12tt \Ba~tl l200ll : \Srwastava et ad l2010l : 
\Zaaarashvili et all |2010| . The QP streaks SI studied here 
exceeded the threshold, but did not show any rapid increase 
in size, brightness, or dramatic change in shape, which means 
that there is no evident signature of the development of the 
kink instability. Moreover, our data process cannot find any 
signatures of the plasma motion along these streaks. Thus, 



these QP streaks could not be the result of the plasma mo- 
tion along twisted helical magnetic field lines. 

A wave train seems to be a more plausible interpretation 
for the QP streaks. A 5-pixel wide slice, labeled 'A' in Fig- 
ure [lid), is placed along the axis of Tube 1 to analyze the 
plasma motion. Figure [21 a) shows the space-time plot gen- 
erated by "seeing" through Slice A in AIA 171 A images. Its 
left vertical axis gives the distance from the start point along 
the slice, and the corresponding height from solar surface is 
marked in the right vertical axis. It could be seen that the 
lower part of the tube oscillated with a period of about 25 
minutes, which was presumably excited by the explosive ac- 
tivity in AR NOAA 11303 at around 09:30 UT (Fig. □»). 
The upward flows mentioned before could be identified as 
stripes with a positive slope in the plot (Fig. EUa)). Linear 
fitting to the stripes suggests that the upward flows moved 
along the tube at a speed of about 32 km s~ . 

The QP streaks SI were located above the distance of 
80 Mm after 10:04 UT in Figure ^ a). Figure ^b) shows 
a zoom-in plot of the blue box region in Figure [21 a) from 
10:03 to 10:18 UT. The temporal evolution of these streaks 
is shown as alternating bright and dark stripes in the space- 
time plot. Obviously, these stripes are steeper than those 
produced by the upward flows. According to the space- 
time plot, the propagation speed of these streaks varied from 
about 70 km s _1 to 120 km s _1 with an average speed of 
about 100 km s . It could also be estimated that the spa- 
tial separation of the streaks is about 5 Mm. Therefore, the 
period is about 50 s if the streaks were caused by a wave 
train. 

Fol lowing the method by \DeForesH [20041 ] and \hiu et aZI 
[201ll ]. Fourier analysis is applied to SI. We generate a three 
dimensional data cube using a series of running difference 
171 A images in the FOV enclosed by the red box in Fig- 
ureQId) from 10:03 to 10:18 UT. The Fourier transformation 
converts the data cube from the space-time domain (x, y, t) 
to the wavenumber- frequency domain (k, cj) (as shown in 
Fig. [2I C ))- A slanted stripe with large powers (painted in 
yellow to red colors) can be clearly seen in Figure EJc). The 
phase speed of the wave can be estimated from the slope of 
the stripe. For that purpose, we equally cut the k-oj plot 
into six pieces along the frequency axis above 10 mHz, and 
identify the point with the maximum power in each piece 
(as marked by the plus signs in Fig. EJc)). The phase speed 
is just the average of the slopes of these points, which is 
about 101 km s _1 with an uncertainty of about 7 km s . 
The integrated power over wave number is presented by the 
black line in Figure EJd), which peaks at about 20 ± 3 mHz, 
corresponding to a period of about 50 ± 7 s. Both results 
are consistent with the direct estimation from the space-time 
plot in Figure HJb) . 

Performing the same analysis to the QP streaks S2 shown 
in Figure [Tie), we obtain a similar result which gives a phase 
speed of about 110 km s _1 and a period of about 45 s. These 
waves were probably excited by the oscillation of the tubes 
or by the upward flows from the tube bottom. 

Recently, some argued that such QP streaks can 
also be explained by h eated QP mass flow s [e.g. , 
De Pontieu and Mclntosfi l2010l: Wclntosh et all |2010l : 



De Pontieu et all l201ll : I Tian efaD. l201l[ |. Those QP mass 
flows may propagate into corona at a speed of order 100 
km s - . To distinguish a wave train from QP mass flows, 
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Figure 2: (a) Space-time plot generated along Slice A from AIA 171 A images with a cadence of 12 s from 09:45 to 10:30 
UT. (b) A zoom-in plot of the region enclosed in the blue box in panel (a) from 10:03 to 10:18 UT. (c) Fourier power 
distribution generated for the 171 A data enclosed by the red box in Fig[lld) from 10:03 to 10:18 UT. (d) Integrated power 
spectrum over wave number of panel (c). (e) Compression ratio in the region between the two green lines in panel (b). (f) 
Average compression ratio over time of panel (e). The red line shows the averaged brightness over time and the green one 
is the fitting line. 



we compare the multi -wavelength observa tions from AIA, 
just like those done bv \Kiddie et all [2012t ] . If the observed 
QP streaks were produced by mass flows, their propagation 
speed would not depend on the temperature, i.e., there is 
no dispersion signature; otherwise, the observed QP streaks 
were caused by a MA wave train. 



The temperature response curves of AIA passbands shown 
in Figure which is generated by the SS W procedure 
' aia_g et .response. pro' provided by AIA team \Lemen et al[ 
|2012| ]. suggest that, except 171 A and 304 A passbands, all 
the other passbands obviously correspond to multiple tem- 
peratures. Thus, 304 A data is used to compare with 171 
A data. It is obvious that there is a similar wave signa- 
ture in 304 A data (Fig. HJa)-(b)), and even more impor- 
tantly, the phase speed is about 70 ± 24 km s" 1 , which is 
much smaller than the phase speed derived from 171 A data. 
However, there is no clear peak frequency in the integrated 
power spectrum (Fig[4jc)); the power roughly stays at the 
same high level from 20 to 40 mHz. This is probably due to 
the relatively low signal-to-noise in 304 A data. Further, a 
comparison of the wave phase between the 171 A and 304 A 
data is given in Figure SJd), in which the time sequences at 
lower, middle and higher positions are presented. The phase 
difference between the two passbands varies with the posi- 
tion changing, which means that the group of QP streaks SI 
is likely a wave train rather than mass flows. The dispersion 
of the MA wave suggests that the plasma tube is likely a 
multi-thread structure, different threads have different tem- 
peratures, and the propagation with different phase speeds 
in different passbands is equivalent to the propagation of 
waves in different threads. 



4 Evidence of slow magneto-acoustic waves 

The QP streaks as viewed in EUV 171 A reflect the fluctu- 
ation of the emission intensity (i.e., brightness), and there- 
fore the fluctuation of the plasma density. It suggests that 
these waves are compressional. In order to obtain the ampli- 
tude or compression ratio of the wave train SI, the parallel- 
ogram defined by the two green dashed lines in Figure Ob), 
where the QP streaks are most clear, is further investigated. 
The background (or unperturbed) brightness in this region 
is a function of distance, which is obtained by simply aver- 
aging the pixel values at each distance, as this region cov- 
ers several wavelengths of the wave. The red line in Fig- 
ure[5Jf) shows the background brightness, /. The amplitude 
of the fluctuation of the emission intensity, dl, is therefore 
obtained by subtracting the background brightness from the 
observed one. In optically thin corona, it is usually true 
that I oc p 2 and therefore the compression ratio, which is 
defined as \dp/p\, could be approxim ately given by \dl /2I\ 
[AschwandeA . \2004 : \Liu et q/.LT201l| . 

Figure Oe) presents the compression ratio in the region 
of interest. Note that the parallelogram is rectified so that, 
along the vertical axis, distance and time variations are cou- 
pled together. Since the green dashed line is selected to 
match the wave speed of 101 km s ~ , the vertical variation 
in Figure [2je) can be roughly interpreted as the variation 
of the wave amplitude for certain phases. The black line in 
Figure [^f) gives the average compression ratio over time, 
which is about 8.4% on average with a standard deviation of 
about 2.0%. No significant damping in amplitude is found 
over the propagation distance of about 33 Mm. Since the 
amplitude of the density fluctuation is relatively small com- 
pared to the background density, it could be treated as a 
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Figure 3: Temperature response curves of AIA passbands. The horizontal dashed lines give the full width at half maximum 
for 171 A and 304 A passbands. 
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Figure 4: (a)-(c) The same to Fig.EJb^d), but for 304 A passband. (d) Time sequences at distance 85, 95 and 105 Mm. 
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linear wave train. In linear theory, of all propagating MHD 
waves, only MA waves can cause density compression. 

The region w e studied i s in a dark cavity, presu mably 
void of plasmas [Z/oml . 1 19961 : \Fuller and Gibson . l2009l | . It is 
hence expected that the magnetic field inside the dark cavity 
is relatively stronger than that outside, so that a pressure 
balance is maintained. Thus the Alfven speed in the dark 
cavity should be approximately equal to or even larger than 
the typical Alf ven speed in qui e t-Sun region, which is about 
220 km s - 1 [Mcintosh et all l201l| . The MA waves we 
investigated here propagated at much slower speeds than 
the Alfven speed, and therefore these waves are slow MA 
waves. 

Since these slow MA waves propagated along the tubes, 
which implies that the wave vector is approximately paral- 
lel to the magnetic field lines, the phase speed of the waves 
should be close to the sound speed. The local sound speed 

is given by c s = *J 2 ^ T , where 7 is the capacity ratio, k 
is the Boltzmann constant, m is the proton mass and T is 
the temperature. The 171 A passband corresponds to the 
temperature of about 0.56 — 1.1 MK, as indicated by the 
horizontal dashed line in Figure [3] The sound speed at this 
temperature range is about 124 — 176 km s _1 . The 304 A 
passband corresponds to the temperature of about 0.07—0.11 
MK, and the sound speed is about 44 — 56 km s . These 
results are close to the phase speeds directly derived from 
the imaging data, but not a perfect match, especially for 
the 171 A data. There could be various explanations for the 
deviation. In a uniform magnetic field, if the wave vector of 
SI was not exactly parallel to the magnetic field, the phase 
speed should be smaller than the local sound speed. Pro- 
jection effect also can result in the underestimation of the 
derived phase speed from the imaging data. Besides, pre- 
vious theoretical analysis showed that, in a magnetic tube, 
slow MA waves propagate at a so called tube speed, which is 
small er than the local sound speed [e.z., \Edwin and Roberts] , 
Il983fl . 

5 Energy flux carried by the waves 

The kinetic energy density of a MA wave can be estimated 
by the formula e = ^pv ph (dp/p) 2 , where p is the plasma 
density and v p h is the phase speed of the wave. The energy 
flux is given by F = ev p h- The analysis of the 171 A data 
in the previous sections has revealed that for wave SI, v p h is 
101±7 km s - 1 and dp/p is 8.4%±2.0%. By further assuming 
the number d ensity to be the typical v alue of 10 s cm -3 in 
dark cavities [Fuller and Gibson . I2009I ]. we infer that the 
energy density carried by wave SI is 3.0 xl0~ 5 ~ 1.2 xlO" 4 
erg cm -3 , or the energy flux is 2.8 x 10 2 ~ 1.1 x 10 3 erg 
cm -2 s _1 . 

Since the emission intensity and therefore the density ob- 
viously decreases with increasing distance (Figure[2jf)), the 
wave energy actually dissipated when it propagated along 
the tube. We applies an exponential equation 7(x) = 
ae c ( x ~ b ' to fit the observed brightness given by the red line 
in Figure [2jf ) . The free parameters a, b and c are found to 
be 434, —0.03 and 80, respectively. The green fitting curve 
suggests that the background brightness weakened by a fac- 
tor of 63% over a distance of about 33 Mm or in about 5 
min. Thus, it is inferred that the average dissipation rate of 



the wave energy is 3.9 x 10 8 ~ 1.6 x 10 7 erg cm 3 s . 

6 Summary 

In this paper, we identified the QP streaks, appearing 
above the quiet-Sun region on 2011 September 25, with the 
slow MA waves. The overall picture of this 'tornado'-like 
event is summarized in Figure [S] There were two flux tubes 
in a dark cavity. Oscillations and upward flows were initi- 
ated in the tubes around 08:00 UT by some perturbations. 
About one hour later, slow MA waves were excited by the 
oscillations and/or the upward flows, and manifested as QP 
streaks in the upper portion of the tubes. The tornado-like 
phenomenon shown in Figure[ljf) was the results of the bulk 
flows along the curved tubes. 

The wave train SI propagated along Tubel around 10:10 
UT is particularly analyzed. Its propagation speed is about 
101 km s _1 in the 171 A bassband or 70 km s _1 in the 
304 A passband, which is comparable with those of the 
slow MA waves found ab ove active regions, 122 km s _1 
[Be Moortel et q/.l. l2002blla| or co ronal holes, 75 - 150 km s 
_1 DeForest and GurmarL I"l998l |. while its period is about 
50 s, which is much shorter than those found above active 
regions, ~ 3 min, or coronal holes, usually 10 — 15 min. 
According to the fluctuation in brightness, the density am- 
plitude of the wave is about 8.4%, also comparable with that 
of the slow MA waves found elsewhere (~ 4% above active 
regions and 10% — 20% above coronal holes, respectively). 
The dissipation rate of the wave energy is estimated from 
3.9 x 10~ 8 to 1.6 x 10~ 7 erg cm -3 s _1 , which is about one- 
order of magnitu de smaller than that o f the r a diative loss of 
the l ocal corona [Withbroe and Noves\ . 1 19771 : [AschwandeA . 
l2006j . 

Our study implies that (1) slow MA waves not only appear 
above active regions or coronal holes, but also exist above 
quiet-Sun regions, and the fact that they were rarely found 
above quiet-Sun region before SDO mission is probably sim- 
ply due to the instrument limitation; and (2) the slow MA 
wave train in this study is not sufficient for local coronal 
heating. It should be noted that the slow MA waves were 
repeatedly observed during the event, and sometimes, two 
or more wave trains appeared simultaneously. Whether such 
slow MA waves ubiquitously exist in the corona and whether 
they do a significant contribution to coronal heating is worth 
to be investigated further. 
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